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\  ABSTRACT 

w 

Since  Rivlin'8  pioneering  work  in  the  non-linear  theory  of  elasticity 
there  have  been  few  new  solutions  which  have  been  exhibited this  work^-we 
establisH^exact  solutions  to  several  boundary  value  problems  in  the  theory  of 
non-linear  elasticity.  A  remarkable  feature  in  all  the  problems  dealt  with  is 
the  possibility  of  an  Infinity  of  exact  solutions. 
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SIGNIFICANCE  AND  EXPLANATION 


Exact  solutions  are  established  for  several  boundary  value  problems  in 
the  non-linear  theory  of  elasticity.  All  the  problems  considered  exhibit  the 
possibility  of  a  multiplicity  of  solutions.  Since  the  domain  involved  in  all 
these  problems  is  infinite,  a  stability  analysis  based  on  energy  criterion  is 
not  readily  available  and  thus  the  stability  of  these  multiple  solutions 
remains  an  open  problem. 
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MEW  EXACT  SOLUTIONS  IN  NON-LINEAR  ELASTICITY 

*  * •*  * 

K.  R.  Rajagopsl  and  A.  S.  Wineman 

1.  INTRODUCTION 

Since  Rivlin’8  pioneering  work  in  the  non-linear  theory  of  elasticity 
(cf.  Green  and  Adkins  [1])  there  have  been  few  new  solutions  which  have  been 
exhibited.  In  this  work  we  establish  solutions  to  several  boundary  value 
problems  in  the  theory  of  non-linear  elasticity.  An  astonishing  feature  in 
all  the  problems  dealt  with  is  the  possibility  of  an  infinite  class  of  exact 
solutions.  The  most  interesting  of  the  problems  studied  is  that  of  the  slab 
in  pure  torsion,  where  in  addition  to  the  classic  torsion  solution  we  exhibit 
the  possibility  of  an  infinity  of  non-symnetric  solutions.  Since  all  the 
problems  considered  involve  infinite  domains,  a  stability  Analysis  based  on 
energy  is  not  readily  available.  We  do  not  study  the  stability  problems  but 
leave  this  to  a  later  date.  The  questions  of  stability  notwithstanding,  the 
exact  solutions  in  themselves  are  so  simple  that  they  are  worth  recording. 

Our  analysis  here  is  an  outgrowth  of  a  recent  study  by  Rajagopal  and 
Wineman  [2]  who  considered  the  deformation  of  a  non-linearly  elastic  material 
which  is  sandwiched  between  two  infinite  parallel  plates,  a  distance  'h1 
apart,  rotated  by  the  same  angular  displacement  fl  about  two  non-coincident 
axes  perpendicular  to  the  plates.  The  form  which  they  assumed  for  the  defor¬ 
mation  correspond  to  a  shear  (which  is  dependent  on  the  co-ordinate  normal  to 
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the  plates)  followed  by  a  rigid  rotation.  They  established  exact  solutions  in 
the  case  of  a  neo-Hookean  and  Mooney-Rivlin  material.  They  also  studied  the 
consequence  of  the  non-linearly  elastic  material  being  characterized  by  a  non- 
convex  stored  energy  function.  In  this  case  they  found  that  it  is  possible 
that  the  solution  possess  discontinuous  deformation  gradients  and  exhibit  the 
phenomena  of  phase  change. 

In  this  work,  we  wish  to  study  associated  but  a  different  class  of 
problems.  All  the  problems  are  concerned  with  the  deformation  of  a  slab  of 
thickness  'h'  whose  other  dimensions  are  infinite,  the  top  and  the  bottom 
surfaces  of  the  sandwich  being  bonded  to  rigid  plates.  Also,  in  all  the 
following  work,  the  non-lxnearly  elastic  material  between  the  plates  is  either 
a  neo-Hookean  or  a  Mooney-Rivlin  material.  First,  it  is  our  aim  to 
investigate  the  consequences  of  rotating  the  top  and  the  bottom  plates  by  a 
constant,  but  differing  amounts,  about  non-coincident  axes  perpendicular  to 
the  plates  (cf.  Figure  1).  For  the  nonlinearly  elastic  materials  in  question, 
we  are  in  a  position  to  establish  an  infinity  of  exact  solution.  He  first 
exhibit  an  exact  solution  in  the  special  case  of  a  neo-Hookean  material.  It 
is  found  that  a  two  parameter  family  of  solutions  is  possible.  He  investigate 
in  detail  the  solution  which  when  the  offset  between  the  axes  'a'  vanishes, 
reduces  to  the  classic  torsion  solution.  However,  this  solution  by  no  means 
is  the  only  solution  to  the  problem  under  consideration.  He  also  compute  the 
tractions  on  the  upper  and  lower  surface  and  find  that  they  are  not  the 
same.  It  is  also  found  that  the  local  contribution  to  the  resultant  moment  is 
different  at  the  upper  and  lower  surfaces  for  fixed  values  of  the  co-ordinates 
in  the  plane  normal  to  the  axes  of  the  plates. 

In  the  case  of  a  Mooney-Rivlin  material,  for  the  kind  of  motion  under 
consideration  the  equations  of  equilibrium  lead  to  a  coupled  system  of  linear 
ordinary  differential  equations  (these  equations  uncouple  in  the  neo-Hookean 
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case).  The  problem  can  be  easily  reduced  to  that  of  solving  a  single  linear 
ordinary  differential  equation  for  a  complex  valued  function.  The  problem  is 
straightforward  and  one  can  determine  the  solution  with  ease.  Once  again  an 
infinite  class  of  solutions  is  possible. 

Next,  we  study  the  problem  of  rotating  the  top  and  the  bottom  plates  by 
different  amounts  about  a  common  axis,  namely  the  torsion  problem.*  Most 
interestingly,  we  once  again  exhibit  an  infinite  one  parameter  family  of 
solutions  for  the  problem  in  question.  The  classic  torsion  solution  belongs 
to  this  class.  While  it  is  the  only  symmetric  solution  to  the  problem,  there 
are  infinitely  many  non-symmetric  solutions  which  are  arbitrarily  close  to  the 
classic  torsion  solution. 

We  also  study  the  problem  of  the  material  being  subject  to  finite  axial 
stretching  followed  by  torsion.  Zn  the  case  of  a  Mooney-Rivlin  material  it 
turns  out  that  the  equations  of  equilibrium  simplify,  after  lengthy  manipula¬ 
tions,  to  a  coupled  system  of  ordinary  differential  equations  which  have  an 
amasingly  simple  solution.  Once  again,  the  problem  exhibits  the  peculiar 
characteristic  of  possessing  an  infinity  of  solutions,  of  course,  when  there 
is  no  stretching,  the  solutions  reduce  to  those  obtained  for  the  pure  torsion 
problem. 

All  the  solutions  in  the  case  of  the  common  axis  problem  and  the  non¬ 
coincident  axes  problem  considered,  except  the  last  class  of  problems,  corre¬ 
spond  to  deformations  in  which  any  plane  parallel  to  the  plates,  remains  in 
the  same  plane,  but  rotates  by  an  amount  ft(z)  which  depends  on  the  coor- 


We  have  received  personal  communication  from  Professor  R.  T.  Shield  [7]  that 
the  above  problem  and  several  other  exact  solutions  are  being  studied 
independently  by  J.  M.  Hill  and  R.  T.  Shield  and  the  work  is  to  be  submitted 
for  publication.  However,  at  this  time  we  have  not  seen  the  manuscript  and 
hence  do  not  know  the  exact  nature  of  their  work. 


dinate  normal  to  the  plates*  That  such  deformations  could  possibly  lead  to 
nonunique  solutions  is  motivated  by  Barker's  investigation  into  the  classical 
Newtonian  fluid  [3] .  That  a  similar  situation  obtains  in  more  general  fluids 
was  established  by  Rajagopal  [4],  [S]  and  Rajagopal  and  Gupta  [6]. 

After  a  brief  review  of  the  necessary  kinematics  in  the  next  section,  we 
derive  the  approximate  form  of  the  equations  of  equilibrium  for  the  problems 
under  consideration  in  section  3.  In  section  4,  we  study  the  problem  of  a 
neo-Hookean  material  sandwiched  between  two  infinite  parallel  plates  which  are 
rotated  about  non-coincident  axes  by  differing  amounts.  The  problem  and 
solution  procedure  in  the  case  of  a  Mooney-Rivlin  material  is  considered  in 
section  5.  In  the  next  section  we  study  the  torsion  problem  in  the  case  of 
the  neo-Hookean  and  the  Mooney-Rivlin  material.  In  the  final  section  we  study 
the  problem  of  the  material  being  subject  to  axial  stretching  followed  by 


torsion 


2.  Kinematics 


For  the  problem  of  rotation  of 
common  axis,  it  would  be  natural  to 

*  ■  tX  -  f (Z)]cos  ft(Z)  - 

y  -  CX  -  f (Z)]sin  ft(Z)  + 

*  ■  z  , 


the  plates  about  differing  axes  and  a 


assume  a  motion  of 

the  form* 

IY  -  g(Z))sin  Q(Z) 

+  flz)  , 

(2.1) 

_ 

_ 

[Y  -  g(Z) }co8  ft(Z) 

+  g(z)  , 

(2.1) 

(2.1) 

where  X,  Y,  z  and  x,  y,  z  represent  the  reference  and  the  current  co¬ 
ordinates  of  the  same  material  point,  respectively.  The  above  motion 
represents  a  deformation  in  which  material  points  which  lie  in  any  plane 
parallel  to  the  plates  continue  to  remain  in  the  plane,  the  plane  rotating 
about  a  point  by  an  amount  fl(z).  The  locus  of  these  centers  of  rotations  is 
in  general  a  curve  in  space  passing  through  the  centers  of  rotation  of  the  top 
and  bottom  plate,  the  locus  being  defined  by 

X»f(Z)  and  Y-g(Z)  .  (2.2),  2 


The  representation  (2»1 >1,2,3  unfortunately  leads  to  certain  singular¬ 
ities  in  the  solution.  It  is  easy  to  see  the  physical  basis  for  this 
difficulty  when  one  considers  that  different  angular  displacements  of  the  top 
and  bottom  plate  could  lead  to  some  plane  parallel  to  these  plates  which  has 
zero  rotation  in  which  case  there  are  difficulties  associated  with  the  notion 
of  the  center  of  rotation.  To  avoid  such  difficulties,  in  this  work  we  choose 
to  express  the  motion  in  the  form 


• 

The  above  form  of  the  motion  is  a  generalization  of  the  motion  assumed  by 
Ra jagopal  and  Wineman  [2]  when  both  the  top  and  bottom  plate  rotate  by  the 
same  aaeunt  ft.  The  form  of  the  motion  assumed  in  [2]  is  the  counterpart  in 
solid  mechanics  of  the  structure  assumed  by  Berker  [3]  and  Rajagopal  [4]  in 
the  case  of  fluids. 


.*  -•  .*  -•  V  _•>  „• 


V  V*  V  .* 


x  -  X  cos  0(Z)  -  Y  sin  Q(Z)  +  f(Z)  , 


y  -  X  sin  Q(Z)  +  Y  cos  0(Z)  +  g(z)  , 


z  ■  Z  . 


(2.3) , 

(2.3) 2 


(2. 3)3 


It  follows  frost  (2.3) ,^2, 3  that  the  defomation  gradient  F  has  the 


following  statrix  representation : 


(£)  -  S 


where 


-xso'  -  yen'  +  f ' 


XCQ'  -  yso'  +  g' 


(2.4) 


C  =  cos  Q(Z)  -  cos  Q(z)  , 


S  =  sin  Q(Z)  ”  sin  fl(z)  , 


(2.5) , 

(2.5) 2 


and  the  prims  denote  differentiation  with  respect  to  the  argument.  One  can 
now  rewrite  the  matrix  (2.4)  in  terms  of  the  current  co-ordinates  x,  y  and 


s  as 


(£)  “  s 


-Q'(y  -  g)  +  f' 
0'(x  -  f)  +  g' 


(2.6) 


For  the  sake  of  simplicity  let  us  introduce  functions  a(y,z)  and  f)(x,z) 


through 


o(y,z)  -  f  -  fl'(y  -  g) 


$(x,z)  -  g'  +  0'(x  -  f) 


(2.7) 


(2.8) 


It  follows  from  (2.6)  -  (2.8)  that  the  left  Cauchy  green  strain  tensor  B  is 
given  by 


.•  V  V  s  «  ’  - 


*  .*  V  -  -  s'  *  ’  -  * 


A'VaV’.  \ 


T  ■  -pi  +  |iB 


3.  Equations  of  equilibrium 

Xn  the  case  of  neo-Hookean  material,  the  equation  of  equilibrium 

div  J  +  pfe  -  0  , 


reduces  to 


3o  A  3 


■  i  h  (*  +  +  If  +  If  “  0  ' 


9 

(3.1) 

:  (.6)  ♦  ||  -  o  . 

(3.2) 

« If  -  H  - 0  • 

(3.3) 

(3.4) 

by  virtue  of  (2.9).  In  equation  (3.1)  div  denotes  the  diverg  e  operator, 
p  the  density  and  |>  the  specific  body  force  field.  In  deriv.  -r  3.2)  - 
(3.4)  we  have  assumed  that  the  body  force  field  is  conservative  and  hence 
derivable  from  a  potential,  i.e.,  £  ”  -grad  $. 

It  follows  from  (2.7)  and  (2.8)  that  (3.2)  -  (3.4)  simplify  to 

-  i  c  If-  0  <  (3-5 


-  i  Jy  <p  *  »♦>  ♦  «•«  ♦  II  • 0  • 


"  u  fit  (p  *  "  0 


It  follows  from  (3.5),  (3.6),  (2.7)  and  (2.8)  that 

32q  w  32g 
dydz  3x3z 

The  equations  (3.9),  (2.7)  and  (2.8)  inq>ly  that 


(3.5) 

(3.6) 


(3.7) 


(3.9) 


and  thus 


-  8"  -  0"  , 


Q(z)  =  ifz  +  , 


where  9  and  are  constants.  The  conditions  that 


fl(-h)  -  Q1  , 


(3.10) 


(3.11) 


...... 


v.'-. 


(3.12) 


fl(h)  =  n2  , 


then  imply  that 


°1  *  °2  ,  ,  “2  - 
2  “a  *  ’  2h 


Next,  equations  (3.5)  -  (3.7)  imply  that 


hr  ♦  !?>  -  0  • 


h  <“'» +  H1  ■ 0  • 

It  then  follows  from  (2.7),  (2.8),  (3.10),  (3.14)  and  (3.15)  that* 


(3.13) 


(3.14) 


(3.15) 


f"  +  ITf  -  0  , 


g"'  +  <*2g'  -  o  . 


(3.16) 


(3.17) 


The  boundary  conditions  which  are  appropriate  for  the  problem  of  rotation 


of  the  top  and  bottom  plates  about  non-coincidents  axes,  are  (cf.  Figure  1) 


f (h)  -  |  sin  fl2,  f(-h)  -  -  |  sin  , 


(3.18) 


g(h)  «  i  ( 1  -  cos  Q2 ),  g(-h)  -  -  j  ( 1  -  cos  ni ) 


(3.19) 


The  above  boundary  conditions  are  however  insufficient  to  determine  the 
solutions  of  equation  (3.16)  and  (3.17).  We  now  proceed  to  obtain  the 
additional  boundary  conditions. 

Integration  of  (3.16)  and  (3.17)  yields 

f"  +  *2f  -  q  , 


(3.20) 


g"  +  f  g  -  r  , 


(3.21) 


where  q  and  r  are  constant.  We  shall  proceed  to  seek  a  solution  in 
which  q  -  r  *  0  which  when  the  offset  between  the  axes  tends  to  zero  reduces 


In  our  analysis  we  shall  restrict  ourselves  to  the  case  when  f  0. 

When  ♦  *  0,  Oj  *  fl2  and  this  case  has  been  studied  in  detail  for  a  much  more 
general  class  of  materials  (materials  with  non-convex  stored  energy  functions) 
and  the  Mooney  material  in  [2]. 


to  the  classic  torsion  problem.  Thus 


f"  +  iff  -  0  ,  (3.22) 

g"  +  *2g  -  0  .  (3.23) 

Equations  (3.22)  and  (3.23)  inqsly  that 

f"(-h)  -  -t2f(-h)  -  a  sin  Q1  ,  (3.24) 

g"(-h)  -  -t2g(-h)  «  a  j-  (1  -  cos  .  (3.25) 

In  the  next  section  we  shall  exhibit  exact  solutions  to  (3.16)  and  (3.17) 
subject  to  the  boundary  condition  (3.18),  (3.19),  (3.24)  and  (3.25). 


Exact  Solution  for  a  Neo-Hookean  Material  in  the  Case  of  Rotation  About 
Non-Coincident  Axes 


For  the  sake  of  convenience  ./id  without  loss  of  generality,  we  shall  set 
»  0  and  “  8q«  this  case  it  is  fairly  straightforward  to  verify 
that  the  solution  to  the  problem  under  consideration  is 


.  a  .  ,0.z+h., 

f(z)  =*  2  sin  — )]  , 


(4.1) 


and 


a  Q0  «0  fl0 

g(z)  ■  j  {[cosec  — - cos  ~  cot  ^p]sin  <(tz 


fl0 

+  sin  —  tan  —  cos  <>z}  , 


(4.2) 


•a  . 

,  .  0 
where  ♦  -  —  . 

It  is  interesting  to  note  that 

fl 

f  (0)  -  |  sin  ~  ,  (4.3) 


and 

n  Q 

g(°)  -  |  sin  Y~  tan  j-  ,  (4.4) 

and  hence  the  functions  f(z)  and  g(z)  do  not  pass  through  the  origin.  It 
is  also  worth  noting  that  when  a  ■  0,  the  above  solutions  (4.1)  and  (4.2) 
reduce  to  the  classic  torsion  solution. 

When  q  and  r  are  non-zero,  one  could  proceed  in  a  similar  manner  to 
establish  exact  solutions.  The  functions  f(z)  and  g(z)  have  the  structure 

f(z)  ”  A  sin  ifz  +  B.  cos  ifz  +  ,  (4.3) 

* 


g(z)  “  A  sin  4>z  +  B  cos  tyz  +  —  ,  (4.6) 

v 

where  the  constants  A^,  A2,  B^,  b2  are  determined  from  the  boundary 
conditions.  Thus,  there  is  an  infinity  of  solutions  possible,  the  solutions 


being  parametrized  by  q  and  r.  The  determination  of  the  solutions  is 


straightforward  and  we  shall  not  provide  the  details  of  the  same. 

Next,  we  determine  the  pressure  field  p.  It  follows  from  equations 
(3.5)  -  (3.7),  (3.20)  and  (3.21)  and  the  definitions  of  a  and  6,  that 


~  U  3x  (?  +  x  +  q  *=  0  , 


(4.7) 


-  -  (p  +  -  *2y  +  r  -  o  , 


(4.8) 


-£fccp  +  *)  -° 


(4.9) 


Thus 


■jj  (p  +  P$)  *  qx  +  ry  -  jjp-  (x2  +  y2)  +  c  . 


(4.10) 


Hence  the  pressure  field  p  is  independent  of  z  if  $  is  independent  of 


We  now  confute  the  tractions  on  the  top  and  bottom  plate  associated  with 
the  above  deformation.  It  follows  from  the  above  solutions  (4.1),  (4.2)  and 


(2.15),  that 


P  I  =  -  ♦yy  + 

“I,  -  th  2 


a  0  .1  y. 

%  2  2  ”  a^  ' 


(4.11) 


-  SiE  (T  t.n 


*  *  +  h 


a  ^2  x  -  1  n0 

11  h  T  +  2  tan  r]  * 


(4.12) 


2  2  Vo 

Thus,  the  resultant  shear  T  ■  (T  +  T  )  *•  on  the  upper  and  lower  surface 

zx  zy 

is  qiven  by 


2  .  .  i, 

T<h)  *  “  i  fI5  *  i  “c  r  ■  + ! tan  r” 2  ■ 


(4.13) 
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(4.14) 


(2  2.  „  n  fl.  y. 

T(-h)  ■  V  *  “  f^T  +  4  sec2  +  <*  tan  - - *)] 


h  2  2  4 

a 


Note  that  the  resultant  shear  traction  on  the  upper  and  lower  plates  differ. 
Also,  when  a  ■  0 


0. 

_  U.  _2 


(4.15) 


c  -  +  h 


(4.16) 


which  corresponds  to  the  "regular  torsion"  problem.  However,  when  q  f  0, 
r  y  0,  the  problems  will  not  reduce  to  the  regular  torsion  problem  (cf. 
section  6). 

The  local  contribution  to  the  resultant  moment  is 


MIS.)  -  <*rEy  -  ,<*„>  lc  ♦  <yT„  +  hTty)l 
* (thT„  -  • 

Thus,  if  Mz  denotes  the  z-component  of  the  local  moment  M, 
M  (+h)  -  (XT  -  yT  )l 

*  *y  «  l+w 


nft  2  2 

+  M  r  f-  -  i  *  i tan  t)]  • 

a 


(4.17) 


(4.18) 


The  first  term  on  the  right  hand  side  is  the  term  corresponding  to  the  usual 
torsion  problem  and  when  a  *  0,  the  expression  reduces  to  the  classic  one. 
It  is  also  interesting  to  note  that 


ME(+  h)  y  Mz(-  h)  when  a  ?  0  . 


(4.19) 


•  .  •  -  *  m  •*  *  *  .  •.  *.*  ’/  *  •  V  V  *  *  * 
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equations 


Introducing  a  complex  valued  function  F(z)  through 


F(c)  -  f(a)  +  ig(z)  , 


(5.11 


we  can  re-write  the  coupled  system  (5.9)  and  (5.10)  as 
F"'  -  itp(j  ~  3$)F'*  -  *2 (V2  -3?)F*  - 
The  appropriate  boundary  conditions  are 


(5.12 


F(h)  -  —  (sin  ♦  id  -  cos  flQ)]  , 


(5.13 


F(-h)  ■  -  —  [sin  8^  +  i(1  -  cos  0^)]  -  0,  since  ■  0  .  (5.14 

As  before,  we  can  determine  the  additional  boundary  conditions  by  integrating 


(5.12).  It  follows  from  (5.12)  that 


F"  -  i*(f  -  3?)F'  -  92 <V2  -  3?)F  -  k  , 


(5.15 


where  k  is  a  oonqplex  constant.  Again,  as  before,  we  can  obtain  an  infinite 
set  of  solutions  parameterized  by  k,  to  the  above  problem. 


The  solution  to  (5.15)  is  of  the  form 

V  *2* 

F(z)  -  Ae  +  Be 


where  Sj  and  s2  are  given  by 


♦2<|  -  3?) 


it  and  »2  -  (1  -  6$) 


(5.16 


(5.17 


The  constants  A  and  B  can  be  obtained  by  using  the  boundary  conditions 


(5.13)  and  (5.14) 


6.  The  Torsion  Problem 

In  this  section  we  consider  the  classic  torsion  problem  and  show  that  an 
additional  infinite  class  of  solutions  is  possible  for  this  problem.  We  shall 
consider  the  problem  of  an  elastic  layer  of  a  neo-Hookean  material  contained 
between  two  infinite  parallel  plates.  The  top  plate  is  rotated  by  an  amount 
Sq,  while  the  bottom  plate  is  held  fixed  (cf.  Figure  2).  For  the  problem  in 
question,  we  once  again  seek  a  motion  of  the  form  (2.3) The  well  known 
torsion  solution  corresponds  to  f(Z)  ■  0  and  g(Z)  «  0,  with  Q(Z)  - 

♦l*  +  V 


Figure  2 

Domain  of  Deformation 


In  the  case  of  a  neo-Hookean  material,  the  equations  of  equilibrium  once 
again  reduce  to 

f"  +  *2f'  -  0  ,  (6.1) 

and 

g'"  +  *2g*  -  0  .  (6.2) 
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■-.y. 


k 


s.y. 

v.v; 


■-  -  -r. 


aru.  a 


The  appropriate  boundary  conditiona  for  the  problem  under  consideration  are 


f  (h)  -  0,  f(-h)  -  0,  g(h)  -  0,  g(-h)  -  0  .  (6.3) 

As  before,  the  above  boundary  conditions  are  not  sufficient  to  determine  the 

solution  to  (6.1)  and  (6.2).  The  additional  conditions  which  can  be  employed 

* 

are 

f(0)  -  A,  g(0 )  -  0  .  (6.4) 


While  the  above  boundary  conditions  can  be  viewed  as  a  kind  of  a  normalization 
of  the  function  f(z),  the  above  conditions  have  a  physical  significance  when 
the  corresponding  functions  f(Z)  and  g(z)  (cf.  equation  (2.2))  are 
considered  as  locus  of  rotations. 

It  is  straightforward  to  show  that  when  th  t  2nw ,  the  solution  to  (6.1) 

and  (6.2)  subject  to  (6.3)  and  (6.4)  is 

m,m\  _  Mcos  6z  -  cos  »h) 

1  '  *  '  14 _ _ J.U  %  9 


g(s)  =  0  . 

When  A  •  0,  we  obtain  the  classic  solution  for  the  torsion  problem.  When 

*  A 

♦h  “  2nv,  one  can  easily  show  that  f(s)  “  A  sin  fz,  g(z)  =  0  where  A  e  R. 

A 

Once  again,  when  A  ■  0,  we  once  again  obtain  the  torsion  solution. 

In  the  case  of  a  Mooney-Rivlin  material,  the  equation  of  equilibrium 
governing  the  torsion  problem  is  (5.15).  The  appropriate  boundary  conditions 
are 


F(h)  -  0,  F(-h)  -  0  and  F(0)  -  A 


(6.6) 


We  can  always  pick  the  axis  in  such  a  manner  that  the  above  is  possible  (cf. 
[3],  [4D. 
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7.  tt»  Problem  of  Finite  Extension  and  Torsion 

In  this  section  we  consider  the  problem  of  a  finite  extension  of  X  of 
the  infinite  sandwich  followed  by  a  twisting  of  the  top  plate  by  an  amount 


Q.  In  this  case,  we  shall  study  a  deformation  of  the  following  kinds 


1 


x  "  — —  [X  cos  Q(XZ)  -  Y  sin  Q(Xz)J  +  f(Xz)  , 
/X 


(7.1) 


1 


y  -  —  (X  sin  a<Xz)  +  y  cos  8(XZ)  , 

/X 


(7.2) 


z  -  XZ  .  (7.3) 

When  X  ■  1,  the  above  deformation  reduces  to  the  expressions  (2.3) , 3 
which  we  have  already  studied  in  detail.  A  simple  confutation  yields  the 
deformation  gradient  £s 


c 

_g 

/X 

✓x 

S_ 

c_ 

/X 

✓x 

0 

0 

✓x(-xsn*  -  yea')  +  f'X 


✓X(XCO'  -  y80')  ♦  g'X 
X 


(7.4) 


where  now  C  “  cos  Q(Xz)  and  S  ■  sin  8(XZ). 


(7,5>1 ,2) 


( 


Z  “ 


V 


in 

terms  of 

C_ 

I§ 

✓x 

/X 

s_ 

C_ 

✓x 

✓x 

0 

0 

(7.6) 


As  before,  we  shall  define 


o(y,z)  =  f'  -  0'(y  -  g)  , 
$(x,z)  a  g»  +  n'(x  -  f)  . 


(7.7) , 

(7.7) 2 


Thus, 
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‘aV-'A 


/  a* 


and  the  left  Cauchy-Green  strain  tensor  £  is  given  by 


and 


Zt  then  follows  that 

IB  -  f  ♦  X2(a2  ♦  62)  ♦  X2  . 

Zt  follows  from  the  constitutive  expression  (2.13)  for  the  Mooney-Rivlin 
Material  and  the  equation  of  equilibrium  (3.1)  that 

“  IS  n  [w(i +  *,(x  +  x2«2,]  *  jg  [“<j ♦ 

+  tv(j  +  J)X2a  ♦  ji(j  -  0)Xa]  -  0  , 

-{*♦!;  ♦  I)X2a$]  ♦  |y  [11(5  ♦  f)(j  +  xV)] 

+  It  +  f)X2$  +  ud  -  P)X5]  -  0  , 


JiS 


« 


”  Pas  3x  +  ®  +  w*2  “ 


+  J7  I»*(j  +  P>x20  +  »<|  -  B>X&1 
-  |j  [y{|  -  p)(^j  +  Xa2  +  Xp2)]  -  0  . 


(7.14) 


It  follows  from  (7.12)  -  (7.14)  and  the  definitions  (7.7),  and  (7.7)2  *or 


a(y,z)  and  P(x,z),  that 

-  |f  +  V(|  +  P)X2P(-0')  +  Wt(j  +  P)X2  +  4  -  P )X]  ff 


0  ,  (7.15) 


“  +  W(j  +  P)X2aft*  +  wl(j  +  P)X2  +  (|  -  P)X]  |f  -  0  ,  (7.16) 


“  If  "  M(5  “  e,x  It  (“2  +  $2)  “  0  • 

It  follows  from  equations  (7.15)  and  (7.16)  that 

32q  _  a2p 

3y3z  3z3x 


(7.17) 


(7.18) 


Thus,  on  using  the  definition  (7.7).,  and  (7.7)2  for  a  and  P,  and  (7.18)  we 


obtain  that 


Thus,  as  before 


«•'  -  0  . 


Q(z)  -  fz  +  • 


(7.19) 


(7.20) 


(tote  that  the  above  expression  for  Q  is  in  terms  of  the  current  co-ordinate. 
Next,  it  follows  from  (7.15)  and  (7.17)  that 

-  11 4  +  f)X2*  +  ji[(~  +  P)X2  +  (I  -  P)X]  ^-2 


+  w(|  -  P)X  yfyj  (a2  +  p2)  -  0  . 

By  virtue  of  (7.7) ^  and  (7.7>2#  (7.21)  can  be  expressed  as 

2 

l“(X2  +  X)  +  P(X2  -  X)]f"  -  t2[(X  -  -  i(2X  +  x2)]f’ 

♦  (f  -  3p)X*g"  -  0  . 


(7.21) 


(7.22) 


V  V1-V-  .-A.- 


-  v  .•  •  -  ■•w  • 


Similarly,  (7.16)  and  (7.17)  inqsly  that 


U(|  +  I)X2t  ||  +  UC<J  +  5)X2  +  -  5)X] 

**  3z 

(7.23) 

*  -  »<J  -  *>»  raj  ♦  «2»  • 


Equations  (7.7)^  and  (7.7)2  now  iroply  that 

tj  (X2  +  X)  +  8(X2  -  X)]g'  *  *  +  ^»2  C^-  X2  -  X  +  ?<X2  +  2X) ]g' 

+  f’tXl-  |  +  3?]  -  0  . 


(7.24) 


Whan  X  *  1,  equations  (7.23)  and  (7.24)  reduce  to  (5.9)  and  (5.10)  for  the 
simple  torsion  problem.  Once  again  we  introduce  a  complex  valued  function 
F(s)  through 


F(z)  «  f(z)  +  ig(s) 

It  then  follows  from  (7.23)  and  (7.24)  that 

2 

[\  (X2  +  X)  +  5(X2  -  X)]F*  * '  -  t2[X  -  -  F(X2  +  2X)]F* 

•  -  i*X(|  -  3$ )F"  -  0  . 


(7.25) 


The  appropriate  boundary  conditions  are 

F(h)  -  0,  F(-h)  -  0  .  (7.26)1<2 

and  F(0)  «  *  .  (7.27) 

In  choosing  the  boundary  conditions  (7.26)j2  and  the  condition  (7.27)  we  are 
basically  choosing  to  interpret  x  ■  f(z)  and  y  *  g(z)  as  determining  the 
locus  of  the  centers  of  rotation.  However,  we  have  indicated  before,  there  is 
no  need  to  make  such  an  interpretation  (indeed  such  an  interpretation  leads  to 
difficulties  if  at  a  plane  z  *  constant,  fl(z)  >0).  One  could  in  fact 
replace  the  condition  (7.27)  by 


F(0)  -  K  =  xQ  +  iyQ 


(7.27)  1 


where  K  is  a  complex  constant  and  in  this  case  one  obtains  a  two  parameter 
family  of  solutions. 

It  is  lengthy  but  straightforward  to  determine  the  solution  to  (7.25) 
subject  to  the  boundary  conditions  (7.26;  .  ^  and  (7.27).  The  solution  is  of 


the  form 


S.z  SjZ 

F(z)  =  Ae  +  Be  +  C  , 


where  Sj  2  are  the  roots  of  the  quadratic  equation 


l\  (X2  +  X)  +  ?(X2  -  X)]s2  -  i*X(|  -  3P)S 


-  4»2[X  -  j-  -  6<X2  +  2X ) ]  «  0  . 


i*X(|  -  36)  +  {-*V(f  -  36 ) 2  +  4 1~  (X2  +  X)  +  6(X2  -  X)] 


(7.28) 


2  1 

[X  -  ~  -  6(X2  +  2X)]*V2 


2  (X2  +  X)  +  B(X2  -  X)]  . 


(7.29) 


Interestingly,  the  above  expression  for  2  simplifies  to 


S,  -  ml  y  - ,) 

‘  X  +  X  +  2B(X  -  X) 


]  -  i* 


(7.30) 


S,  •  ml  1  -  ”'V  *’  1  • 


X2  +  X  +  26 (X2  -  X) 


(7.31) 


It  then  follows  that 


f(Xz)  m  —  tsin  $hXM  sin  <|»Xz  -  sin  ♦Xh  cos  i|>MXZ  +  A] 
A 


(7.32) 


g(Xz)  ■  —  tsin  <J>hXM  sin  i|»Xz  -  sin  \|>Xh  sin  ^MXz]  , 
A 


where 


(7.33) 


A  -  sin  »hX(1  -  M)  ,  M  -  ]lhH\ fr-f  .  (7.34) 

and 

A  -  4  sin  *hX  (-  “ sin  ^  sin  .  (7.35) 

We  conclude  this  work  by  studying  the  problem  of  the  infinite  sandwich 
being  extended  by  a  finite  amount  and  then  being  twisted  by  0^  and  fl2 
about  non-coincident  axes*  For  this  problem,  we  once  again  seek  a  deformation 
of  form  (7.1)  -  (7.3).  Thus  we  obtain  the  same  equilibrium  equation  (7.25). 
The  only  difference  is  in  the  boundary  condition  and  for  the  above  problem, 
the  boundary  conditions  are  (for  the  sake  of  convenience  we  shall  set  0^=0 
and  Q2  -  fl0) 

f  (Xh)  -  —  sin  fl,  -  sin  Q.  ,  (7.36) 

2/X  2  2/X  0 


and 


Thus 


f  (-Xh)  -  -  “  sin  fl  «  0  , 

2/\  1 


g(Xh)  *  (1  -  cos  fl. )  »  (1  -  cos  R.)  , 

2/X  2  2/X  ° 


g(-Xh)  ■  — ( 1  -  cos  fl.)  -  0  . 

2A  1 


F(Xh)  ■  tain  fl.  +  id  -  cos  fl.)]  , 

a  /r  "  u 


F(-Xh)  -  0  . 


(7.37) 

(7.38) 

(7.39) 

(7.40) 

(7.41) 


Since  fl.  ■  0  and  fl.  ■  fl..  it  follows  that 


Note  that  the  above  yield  a  one  parameter  family  of  solutions  in  terms  of  the 


constant  A3,  and  this  is  to  be  expected  as  we  have  not  imposed  an  additional 
condition  as  before,  e.g.  the  location  where  the  locus  of  the  centers  of 
rotations  cut  the  z  -  0  plane. 
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